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experiments currently under construction. We investigate synergies between the 
NOvA long-baseline accelerator experiment with atmospheric neutrino data from 
the India-based Neutrino Observatory (INO). We identify the requirements on en- 
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1 Introduction 

Huge progress has been achieved in the study of neutrino oscillations [1-5] and a rough 
picture of the structure of three-flavour lepton mixing has been obtained, with two large 
mixing angles {612 and ^23) and two neutrino mass-squared differences separated roughly 
by a factor 30: Amgi ~ 7.6 ■ 10~^ eV^ and |Am|]^| ~ 2.4 ■ 10"^ eV^, see [6,7] for a recent 
global fit. The sign of Am^i is determined by the matter effect [8-10] inside the sun being 
responsible for the flavour transition of solar neutrinos.^ In contrast the sign of Ato^^ is 
not known and present data cannot distinguish between the so-called normal or inverted 
neutrino mass ordering, with Am'^i > or Am|^ < 0, respectively. The determination 
of the sign of A?n|^ is one of the most important goals of the future neutrino oscillation 
program. The type of the neutrino mass ordering provides crucial information on the 
flavour structure in the lepton sector. Furthermore, the unknown sign of Amg^ introduces 
a two-fold ambiguity in the determination of the parameters ^13 and 6 by long-baseline 
experiments [11, 12], and might severely affect the search for CP violation in neutrino 
oscillations. 

The most promising way to determine the neutrino mass ordering is to search for 
matter effects [8-10] in oscillations driven by Amg^^. This requires the participation of 
the electron neutrino in such oscillations, which is suppressed by the mixing angle ^13. 
Therefore, the size of ^13 strongly affects the possibility to determine the sign of Aml^- 
Recently the Daya Bay reactor neutrino experiment established a non-zero value of 613 at 
the 5cr level [13], 

sin^ 2^13 = 0.092 ± 0.016(stat) ± 0.005(syst) . (1) 

A similar result was later found also by the RENO experiment [14] 

sin^ 2^13 = 0.113 ± 0.013(stat) ± 0.019(syst) . (2) 

These results confirm previous hints for a relatively large value of 613 from the T2K [15] 
and Double Chooz [16] experiments, see also [7]. Given these exciting developments, the 
question arises whether there is a chance to determine the neutrino mass ordering with 
currently running or planned experiments. This possibility will have an important impact 
on the planning and design of a subsequent generation of oscillation experiments towards 
the search for leptonic CP violation. 

The main goal of the current generation of accelerator experiments (T2K [17] and 
NOvA [18]) as well as reactor experiments (Double-Chooz [19], RENO [20], and Daya 
Bay [21]) is the determination of ^13, see [22] for a recent review. From those experiments 
only NOvA may have sensitivity to the sign of Arn^^. However, the study performed 
in [23] shows that, combining expected data from NOvA with the ones from all the other 



^The convention-independent statement is that the neutrino mass state which contains dominantly the 
electron neutrino (denoted by vi per convention) has to be the hghter of the two mass states responsible 
for the flavour transitions observed for solar neutrinos and reactor anti-neutrinos in the KamLAND 
experiment. 



mentioned experiments, only a poor sensitivity to the neutrino mass ordering will be 
obtained: for sin^ 2^13 ?a 0.09 the sign of Am|^ can be determined at 90% CL in 2019 only 
for about 45% of all possible values of the CP phase 6 and there is negligible sensitivity 
at higher CL. Even for a fully optimised run time schedule as well as optimistic upgrade 
assumptions for T2K and NOvA, the global data from those experiments extrapolated 
until 2025 will provide sensitivity at 3a for only about 35% of all 6 values. For similar 
studies see also [24-26]. 

Motivated by this situation we explore here the possibility to use data from atmo- 
spheric neutrinos in order to determine the mass ordering. In particular, in the India- 
based Neutrino Observatory (INO) [27] a magnetized iron calorimeter will be built for the 
observation of charge separated muons induced from atmospheric neutrinos, with a 50 kt 
detector expected to start data taking in 2017 [28]. A comparable small sample of such 
atmospheric neutrino events has been observed by the MINOS detector [29,30]. For large 
enough values of 6*13 matter effects will induce characteristic signatures for atmospheric 
muon neutrinos as a function of the neutrino energy and zenith angle, different for neu- 
trinos and anti-neutrinos, depending on the mass ordering, due to three-flavour matter 
effects [31-34], see also [35]. The sensitivity of magnetized iron detectors to those effects 
has been studied in [36-44]. Also non-magnetized atmospheric neutrino detectors may 
provide sensitivity to the neutrino mass ordering, see e.g., [43,45-50]. However, the lack 
of event-by-event discrimination of neutrino and anti-neutrino induced events leads to a 
dilution of the relevant signatures and therefore huge detectors will be required [51,52]. 

In the following we will investigate the impact of the INO atmospheric neutrino data on 
the sensitivity to the neutrino mass ordering in a global context by combining simulated 
data from INO, NOvA, and T2K. The remainder of this work is organized as follows. 
In section 2 we give details of our simulation of atmospheric and accelerator data, as 
well as on the statistical analysis. The main results are presented in section 3, where we 
discuss the combined sensitivity of atmospheric and accelerator data as a function of time 
for various assumptions on the experimental configuration achieved in INO. In section 4 
we discuss the effect of neglecting the solar mass-squared difference Amli in the INO 
analysis. We conclude in section 5. 

2 Simulation details 

2.1 Atmospheric neutrinos in INO 

For the simulation of atmospheric data in INO we follow closely [40], where technical 
details for the calculation of the event rates, the used cross section and neutrino fluxes, 
as well as the statistical analysis are given. Here we summarize our main assumptions. 
We assume a muon threshold of 2 GeV and assume that muon charge identification is 
perfect with an efficiency of 85% above that threshold. As stressed in [39,40] the energy 
and direction reconstruction resolutions are crucial parameters for the sensitivity to the 



mass ordering. The ability to reconstruct neutrino energy and direction depends on the 
corresponding resolutions for the muon, the mean angle between muon and neutrino, 
and the energy and momentum reconstruction for the hadronic shower. In the absence of 
detailed Monte Carlo simulations we assume that neutrino energy and direction resolutions 
are Gaussian with widths aE and erg, respectively, assuming two representative sets of 
values, corresponding to a "low" or "high" resolution configuration: 

(Tij/E, = 0.15, ae = lh° (low) 

(Te/E, = 0.10, ae = 10° (high) 

We take those resolutions to be independent of energy or zenith angle. In a real experiment 
neither the resolutions nor charge ID efficiencies will be constant, and one may expect 
different data samples with diverse reconstruction qualities. Unfortunately detailed re- 
construction capabilities of the INO detector are currently not available. Therefore we 
make the simplified assumptions stated above. Our results should be considered as a 
representative estimate assuming that the adopted values can be achieved in average for 
the majority of the events. With a constant efficiency of 85% and a muon threshold of 
2 GeV we find 242 (223) /i-like events per 50 kt yr exposure for the high (low) resolution 
assuming no oscillations (sum of neutrino and anti-neutrino events). 

As shown in [40] e-like events provide additional sensitivity to the mass ordering. 
However, electrons and positrons induce electromagnetic showers in an INO-like detector 
and the reconstruction and charge separation of such events is difficult. Currently it is 
not foreseen to consider e-like events in INO. Therefore we do not include them in our 
analysis. 

We divide the simulated data into 20 bins in reconstructed neutrino energy from 2 GeV 
to 10 GeV, as well as 20 bins in reconstructed zenith angle from cos ^ = — 1 to cos 9 = —0.1. 
We then fit the two-dimensional event distribution in the 20 x 20 bins by using the 
appropriate x^-definition for Poisson distributed data, while also including the following 
systematic uncertainties: we assume a 20% uncertainty on the over-all normalization of 
events, and 5% on each of the neutrino/anti-neutrino event ratio, the z/^ to u^ flux ratio, 
the zenith-angle dependence, and on the energy dependence of the fluxes, see [40] for 
details. 

In the simulation of atmospheric neutrino data we set Am'^i = 0. This is a rea- 
sonable approximation for neutrino energies above 2 GeV. We estimate the accuracy of 
this approximation in section 4. Under this assumption, atmospheric neutrino data de- 
pend on the three parameters Am'^Q, 623,613. For the NOvA and T2K simulations we 
use full three-flavour oscillation probabilities including Atti^i effects. In order to combine 
such simulations with INO we use for the atmospheric analysis an effective mass-squared 
difference, which is related to Aml^ by [53] (see also [54]) 

Am^g = Amg^ — (cos^ 6^12 — cos 6 sin 6^13 sin 26^12 tan 623) Aml^ . (4) 

This is particularly important when investigating the sensitivity to the mass ordering. 



since to very good approximation the degenerate solutions occur at ±A?Ti^jj, which cor- 
responds to shghtly differing values of |A?7i|^| according to the above equation. 

2.2 NOvA, T2K, and reactor experiments 

For the simulation of the long-baseline accelerator experiments NOvA and T2K we use 
the GLoBES software [55,56] and follow closely the analysis of [23], where details on the 
assumed experimental parameters can be found. For T2K, we always assume the final 
exposure corresponding to 5 years running with a beam power of 0.75 MW in neutrinos 
only. Assuming a realistic beam power evolution, such an exposure will be obtained 
around 2018. For NOvA, the nominal exposure is 3 years for neutrinos and anti-neutrinos 
each, with 6- 10^" POT per year and a 15 kt detector. When considering the time evolution 
of the sensitivity we assume the run plan for alternating neutrino and anti-neutrino beams 
as well as the time evolution of the detector mass as shown in Fig. 11 of [22]. This should 
be considered as one particular representative example, and in this scenario the nominal 
exposure would be reached around 2019 for neutrinos and 2020 for anti-neutrinos. We 
assume that the experiment is terminated after this date. 

The main sensitivity to the mass ordering comes from NOvA due to the somewhat 
larger matter effect, because of the longer baseline of 810 km. The matter effect is 
significantly smaller in T2K (295 km baseline) and therefore T2K has no sensitivity to 
the mass ordering. However, T2K may contribute indirectly by providing additional 
information on 613 and the phase 6. Therefore, we always consider the combination of 
both beams in the following. 

We include the determination of sin^ 2^13 from reactor experiments as a simple Gaus- 
sian prior centered at 0.09 with the current Daya Bay error of 0.017 at la (systematic and 
statistical errors combined). While improved data from the running reactor experiments 
will tighten these bounds, the actual sensitivity of INO and NOvA to the mass ordering 
does not depend crucially on the exact error. 

2.3 Parameter values and implementation 

For our simulations, we have fixed the values of the true neutrino oscillation parameters 

to 

Amis = 2-4 ■ 10-3 eV^ Am^^ = 7.8 ■ 10"^ eV^ 

^23 = 45°, 012 = 33°, sin^ 2^13 = 0.09, ^ ^ 

unless stated otherwise. When treating the accelerator experiments, the sensitivity to 
the mass ordering depends significantly on the true values of the CP phase 6. We fix 
Am^i and Ou in our fits, since the impact of these parameters is expected to be minimal. 
Other parameters are marginalized over, and we impose the following priors {la Gaussian 
errors): a{sm'^ 2613) = 0.017 (the uncertainty from Daya Bay [13]), a{Aml-^^) = 0.5|A?Ti^3| 
(this is a very weak prior with the only purpose to guide the minimization algorithm), and 
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Figure 1: The dependence of the marginahzed atmospheric x^ as a function of Am|j^ in the different 
mass orderings assuming jAm^jf | ~ 2.4 • 10"'^ cV^. The dotted (sohd) curves correspond to sin^ 26'i3 = 
(0.09) and the black (red) curves to a detector mass of 50 kt (100 kt). The running time has been assumed 
to be 10 years with the high resolution scenario. 



(^{(^23) = 0.086'23. The prior on 6^23 corresponds to cr(sin 6*23) ~ 0.063 at ^23 = 45°, which 
is approximately the accuracy from current data [6]. Once T2K and NOvA are included in 
the fit, they will provide a more accurate determination of ^23 than this prior. By default 
we assume that the true mass ordering is normal and test the sensitivity to exclude the 
inverted ordering. With the above mentioned priors the results are very similar if the 
true ordering is inverted, as we will show explicitly towards the end of section 3. For 
each simulated value of the true parameters, we find the minimum value of the x^ iii the 
opposite mass ordering. When reporting sensitivity in terms of standard deviations we 
use one degree of freedom to evaluate the x^- Hence the number of a with which the 
wrong mass ordering can be excluded is given by the square-root of the Ax^ between the 
two signs of Am^^ . 

In Fig. 1, we show the behavior of the marginalized x^ ^^ a function of Arrig^ for the 
case of atmospherics only for sin^ 26'i3 = and 0.09, respectively. As can be seen from 
the figure, the mass orderings are indistinguishable when sin^ 26*13 = 0, while there is a 
significant difference for sin^2^i3 = 0.09. Also note the shifts in the best fit of Am'^i as 
compared to the simulated Amlg, which are well in agreement with Eq. 4. 

Let us stress that the sensitivity of atmospheric data on the mass ordering also depends 
on the value of ^23, with typically improved (weaker) sensitivity for ^23 > 45° (< 45°). 
This effect has been investigated in [40]. Here we always assume a true value ^23 = 45° 
in order represent the "average" sensitivity and briefly show this behavior in the end of 
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Figure 2: The sensitivity of an INO-like detector only as evolved in time assuming start of data taking 
in the beginning of 2017 (ticks in the figure correspond to the beginning of the year). We show the 
number of standard deviations with which the wrong mass ordering can be excluded. Black (red) curves 
correspond to a detector mass of 50 kt (100 kt) and dashed (solid) curves correspond to the low (high) 
resolution scenario (see text). We assume true values sin^ 20^3 = 0.09, 6'23 — 7r/4, and impose external 
priors at la of a{sm^ 26*13) = 0.017 and ct(6'23) = 0.0=^'^ 
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the next section. 



3 Results 



Let us first investigate the sensitivity of atmospheric data alone. In Fig. 2, we show 
how the sensitivity of INO only to excluding the inverted ordering would evolve in time, 
depending on implemented scenario. We would like to remind that to a fairly good 
approximation this sensitivity is independent of S and will be slightly worse if the effects 
of this parameter are included (see Sec. 4). From this figure, we can deduce that the 
sensitivity of atmospheric data would be strongly dependent on how well the atmospheric 
study can be performed. In particular, in comparing the best and worse cases, the high 
resolution 100 kt scenario would reach a 90% CL sensitivity after slightly less than two 
years of running, while the low resolution 50 kt scenario would not accomplish this within 
the assumed 10 year lifetime of the experiment, after which the high resolution 100 kt 
scenario has reached a sensitivity close to 3cr. Note that the high resolution 50 kt scenario 
outperforms the low resolution 100 kt scenario, implying that reaching for an increase in 
resolution may be preferable to increasing the detector mass [39,40], depending on the 
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Figure 3: The minimum and maximum sensitivities (depending on the true value of 5) of atmospherics 
combined with NOvA and T2K. We show the number of standard deviations with which the wrong 
mass ordering can be excluded. The left (right) panel corresponds to a detector mass of 50 kt (100 kt) 
and dashed (solid) curves correspond to the low (high) resolution scenario. The shaded area is the 
corresponding result for NOvA and T2K only. The true value of sin'^ 26*13 has been assumed to be 0.09. 



cost and technical feasibility of doing so. 

The sensitivity to exclude the inverted mass ordering by combining information from 
INO together with NOvA and T2K is presented in Fig. 3. Since the sensitivity depends 
on the true value of the CP violating phase 5, we show both the maximum and minimum 
sensitivities, corresponding to the most and least favourable values of 5, respectively. For 
comparison, we have also included the sensitivity of the accelerator based experiments 
only. The kink in the least favourable curve is due to the onset of anti-neutrino running 
in NOvA, which in our example time line happens simultaneously with the onset of INO 
data taking in 2017. From this figure, we see that in the low resolution scenarios the 
sensitivity to the mass ordering is mildly improved compared to that of the accelerator 
experiments only, increasing by roughly 0.5o" — Icr after the full 10 years of running, 
depending on the scenario and the value of 6. For these scenarios, the atmospheric data 
only slightly adds information to the accelerator data during the running time of NOvA, 
resulting in a given sensitivity to be reached marginally ahead of when it would be reached 
by accelerator experiments alone. 

For the higher resolution scenarios, the rise in the precision is significantly faster, in 
particular during the time when atmospheric and accelerator experiments are running in 
parallel. This results in a significant improvement in the sensitivities in a relatively short 
time span. It should be noted that this effect is present both for the values of 6 to which 
accelerator experiments are most sensitive, as well as those values to which they are least 
sensitive, leading to an overall increase in the sensitivity irrespective of the value of 6. 
The statistical level where the synergy between atmospheric and accelerator experiments 
is most apparent is around the 2a region for the least favourable values of 5, where the 
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Figure 4: The time evolution of the fraction of values of the CP violating phase 5 for which the 
combination of INO, NOvA, and T2K would be sensitive to the mass ordering at 2a. Black (red) curves 
correspond to a detector mass of 50 kt (100 kt) and dashed (solid) curves correspond to the low (high) 
resolution scenario. The shaded area is the corresponding result for NOvA and T2K only. The true value 
of sin^ 2013 has been assumed to be 0.09. 

final sensitivity is impacted by both types of experiment. Most noticeably, NOvA and 
T2K never reach this sensitivity by themselves, while the high resolution atmospheric 
data by itself would reach it several years later than the combined. Most strikingly, the 
combined sensitivity in the most optimistic high resolution 100 kt scenario would reach 
the 2a within 2 years of the start of atmospheric data taking, still within the lifetime of 
NOvA. 
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Figure 5: Same as Fig. 4 but for sensitivity at the 90% (3cr) CL for the left (right) panel. 
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Figure 6: The mininiuin and maximum sensitivities (depending on the true value of S) of atmospherics 
combined with NOvA and T2K assuming that the true mass ordering is normal (inverted) for dashed- 
black (solid-blue). We assume a detector mass of 50 kt, low resolution, and a true value sin^ 2^13 = 0.09. 

In Fig. 4, we show the time evolution of the 2a CP fraction, the fraction of the values 
of true S (assuming a flat distribution) for which the sensitivity to the mass ordering 
is 2a or better. Note that this information is complementary to that presented in the 
previous figure, as it also contains information about the sensitivity dependence on 6, 
not just the maximum and minimum sensitivity. Comparing these two figures, we can 
deduce that unfortunately there are more values of 6 which are close to the minimal 
sensitivity than to the maximal one. In particular, for the worst case of low resolution 
and a 50 kt detector, the 2a CP fraction remains around 0.5 by 2027, even though the 
minimal sensitivity at this point is around 1.9cr. This is further illustrated by the rate at 
which the CP fraction grows after reaching 0.5 for the other scenarios. For completeness, 
we also show the 90% and 3a CP fractions in Fig. 5. We deduce from these figures that 
it would be relatively easy to obtain a 90% CL hint for rejecting the inverted ordering, 
while a 3a evidence will be significantly more challenging. However, in both cases, the 
interplay between atmospheric and accelerator data will allow for establishing these at an 
earlier time, regardless of the true value of 6 or the confidence level in question. 

So far we always assumed that the true mass ordering is normal. In Fig. 6 we compare 
the difference in sensitivity between a true normal or inverted mass ordering, confirming 
that the sensitivity is very similar. The features induced by the matter resonance will 
appear in the data either for neutrinos or anti-neutrinos, depending on whether the true 
ordering is normal or inverted, respectively. While statistics will be larger for neutrinos 
than for anti-neutrinos, the important quantity for the sensitivity is the difference in 
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Figure 7: The minimum and maximum sensitivities (depending on the true value of 6) of atmospherics 
combined with NOvA and T2K as a function of the true value of ^13 (left) and ^23 (right), respectively. 
We show the number of standard deviations with which the wrong mass ordering can be excluded. We 
have assumed the total final exposures for NOvA and T2K, and an exposure of 500 kt yr for INO. 
Dashed (solid) curves correspond to the low (high) resolution scenario. The gray-shaded area is the 
corresponding result for NOvA and T2K only. The blue-shaded areas indicate the current 1, 2, 3a 
regions of the parameters (from dark to light shading) . 

event numbers expected for the different mass orderings in each of the two samples, and 
this difference is largely independent of the true ordering. This statement is true for 
approximately fixed oscillation parameters, especially 6*13. Our simulation shows that 
the accuracy on ^13 and other parameters as provided by current Daya Bay data as well 
as the simulated data from T2K and NOvA, is sufficient for the above argument to be 
valid, in agreement with [40]. For Fig. 6 we have assumed a 50 kt detector mass and low 
resolution for INO. Corresponding results for all other configurations considered in this 
work are very similar to the one shown in Fig. 6. 

Despite the strong evidence for a non-zero (and relatively large) value for 613 the 
current uncertainty from Daya Bay [13] still allows for a significant spread in sin^ 26*13, 
see Eq. 1. In Fig. 7 we show the dependence of the sensitivity to the mass ordering as 
a function of the true value of sin^ 26*13, assuming the final exposure for all experiments. 
The blue shading indicates the 1, 2, and 3a lower bounds on sin^ 2^13 according to Eq. 1. 
We observe that if the true value of 6*13 happens to be close to the current 3a lower bound, 
the sensitivity to the mass ordering will be noticeably weaker than at the current best fit 
point, typically allowing for an exclusion of the wrong hierarchy with about one standard 
deviation smaller significance. Correspondingly better sensitivities can be achieved for ^13 
larger than the current best fit point (not shown). In a similar fashion, we also show the 
dependence of the sensitivity on the mixing angle 6*23. As shown in [40], the sensitivity 
becomes better for values of 623 above ir/A and worse for smaller values, see e.g., Eq. 32 
of [40]. It is also notable that for smaller values the atmospherics help only marginally. 
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Figure 8: A^^ of the wrong mass ordering for sin^ 20i3 = 0.09 including a finite Arn^i in the INO 
simulation as a function of the CP phase S. Data is simulated for normal mass ordering and 6 = 0. Those 
data are fitted assuming the inverted mass ordering varying 6. When changing the sign of Arri^i we 
keep lAwgfjl constant, see Eq. 4, and all other oscillation parameters are fixed. For INO only statistical 
errors are taken into account, corresponding to an exposure of 500 kt yr, and we show results assuming 
high and low resolutions. The horizontal thin lines correspond to the Ax^ for Arn^i = 0. For NOvA we 
assume 3 yr nominal exposure for neutrinos as well as anti-neutrinos. 

while for larger values the additional information from atmospherics have a larger impact. 

4 On the size of Amli effects in INO 

For our simulations of INO data we neglect effects of Am'^i. Oscillations of atmospheric 
neutrinos due to the solar mass-squared difference are important for neutrino energies 
below 1 GeV, and since we always impose a threshold of 2 GeV Am^i effects are expected 
to be small. The approximation Arn^i = greatly simplifies the numerical calculations 
and makes the detailed analysis including parameter correlations, detector resolutions, as 
well as systematical uncertainties feasible. In this section we check the accuracy of this 
approximation by considering a simplified analysis taking into account full three-flavour 
oscillation probabilities. 

We consider only statistical errors corresponding to a 500 kt yr INO exposure and 
neglect systematical uncertainties in the fit. We simulate "data" assuming normal mass 
ordering, sin^ 2^13 = 0.09, and 6 = 0. Those "data" are fitted with inverted mass ordering, 
where we keep |Amgg| constant according to Eq. 4. All other oscillation parameters are 
fixed at the values assumed in generating the data, except for the CP phase 6. The Ax^ is 
shown in fig. 8 for low and high resolutions in INO. We see that x^ changes by about one 
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unit as a function of S. In the figure we also show the Ax^ value obtained by assuming 
Amli = 0. We conclude that full three-flavour effects do not lead to additional sensitivity 
to the mass ordering, but will diminish the sensitivity by less than one unit in x^ when 
marginalized over 6. This is true for both assumptions on the resolution. In this sense 
our INO sensitivities are slightly optimistic. One can expect that those results hold also 
in the presence of systematical effects. 

The various assumptions we had to make at the current stage on reconstruction ca- 
pabilities and efficiencies introduce uncertainties on the sensitivity which are most likely 
larger than the one from neglecting the solar mass-splitting, justifying our approxima- 
tion. However, we stress that once reliable detector properties become available it will 
be necessary to include full three flavour effects in order to obtain accurate sensitivity 
predictions. 

In Fig. 8 we show also the 6 dependence of the x^ from NOvA data for the same 
type of analysis. We see that for NOvA three-flavour effects are essential and there is a 
significant dependence on the CP phase 6."^ Comparing the variation of INO and NOvA 
we conclude that there is very little complementarity with respect to 6 dependent effects. 
The 6 best fit point will be completely dominated by NOvA, and hence the small decrease 
in sensitivity for INO might be somewhat lifted once combined with NOvA. This results 
in x^ values very similar to the one obtained in the approximation Attt.^]^ = 0. 

5 Conclusions 

The recent discovery of a relatively large value of 613 by the Daya Bay reactor experiment 
opens exciting possibilities for the future neutrino oscillation program. In this paper we 
have focused on the determination of the neutrino mass ordering being normal, Am'^i > 0, 
or inverted, Arrigj^ < 0. Currently planned long-baseline accelerator experiments, in 
particular NOvA and T2K, only have a poor sensitivity to the mass ordering even in case 
of large ^13. On the other hand, a large ^13 provides potentially interesting opportunities 
for the INO atmospheric neutrino experiment, which is planning to start data taking in 
2017. We have investigated the combined sensitivity of these experiments for different 
assumptions on the size and event reconstruction capabilities of INO. 

Assuming "low" resolutions (15% and 15° reconstruction accuracy in neutrino energy 
and direction, respectively) and a 50 kt detector (fiducial) only a poor global sensitivity 
is obtained and there is only marginal improvement from combining INO with NOvA and 
T2K data. Significant synergy and improved sensitivity is obtained for "high" resolution 
(10% and 10°) and/or doubling the INO detector mass. We find that improving the 
resolution is more effective than increasing the detector mass. For high resolution we find 
that a 2a determination of the mass ordering is possible irrespective of the CP phase 6 



^Thc NOvA curve in Fig. 8 depends strongly on the true value of S, whereas the result for INO depends 
only weakly on the assumed value. 
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in 2021 (2019) with a 50 kt (100 kt) detector. The high resolution 100 kt detector even 
allows a determination at 3a around 2025. These conclusions hold for sin^ 2^13 = 0.09, 
close to the current best fit point. The sensitivity of the considered experiments still 
depends crucially on the actual true value of ^13 within the currently allowed 3a region, 
and to a lesser extent also on 623. 

In our analysis of simulated INO data we assumed a constant efficiency of 85%, and 
energy and direction resolutions to be independent of energy and direction. These are 
simplifying assumptions and once a detailed detector simulation becomes available a more 
realistic analysis should be carried out. We have also discussed the impact of effects related 
to the solar-mass splitting Am^i and the CP phase 6 for the atmospheric neutrino data 
and have shown that their impact on the sensitivity to the mass ordering is small. 

In conclusion, the by now established large value of ^13 opens the possibility to deter- 
mine the neutrino mass ordering within a time frame of about ten years with experiments 
currently under construction. Atmospheric neutrino data from INO may be crucial in or- 
der to achieve this goal and we believe that it is important to include such synergies in the 
global context towards future neutrino oscillation facilities. Our study suggest, however, 
that in order to achieve a relevant sensitivity some improvements of the INO detector 
(either in event reconstruction capabilities, detector mass, or both) seem to be necessary. 
Complementary information could be provided by e-like events, and the sensitivity of INO 
could be potentially increased significantly if the reconstruction of charge-separated e-like 
events was possible [40]. 

Acknowledgement. We thank Tarak Thakore, Nita Sinha, and Enrique Fernandez- 
Martinez for useful discussions and Anselmo Meregaglia for spotting a typo in Eq. 4 in 
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